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1. Introduction 



An ideal charged boson system is of interest because of the simplicity in 
its formulation and yet the complexity of its manifestations. The astonishingly 
complicated behavior of this idealized system may provide some insight to the 
still not fully understood properties of high T c superconductivity. As is well 
known, R. Schafroth[l] first studied the superconductivity of this model fifty 
years ago. In this classic paper he concluded that at zero temperature T = 
and in an external constant magnetic field H , there is a critical field 



with p denoting the overall number density of the charged bosons and m, e 
their mass and electric charge respectively; the system is in the super phase 
when H < H c , and in the normal phase when H > H c . Due to an oversight, 
Schafroth neglected the exchange part of the electrostatic energy, which inval- 
idates his conclusion as was pointed out in a 1990 paper [2] by Friedberg, Lee 
and Ren (FLR). This oversight when corrected makes the ideal charged boson 
model even more interesting. Some aspects of this simple model are still not 
well understood. 

In what follows we first review the Schafroth solution and then the FLR 
corrections. Our discussions are confined only to T = 0. 



Let 0(r) be the charged boson field operator and <^(r) its hermitian con- 
jugate, with their equal-time commutator given by 



These bosons are non-relativistic, enclosed in a large cubic volume £1 = L 
and with an external constant background charge density —ep ext so that the 
integral of the total charge density 



(H c )sch = ep/2m 



(1.1) 



2. Hamiltonian and Schafroth Solution 



Wr), 0t(r')]=<5 3 (r-r'). 



(2.1) 




(2.2) 



is zero. The Coulomb energy operator is given by 




(2.3) 
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where : : denotes the normal product in Wick's notation[3] so as to exclude 
the Coulomb self-energy. 

Expand the field operator 0(r) in terms of a complete orthonormal set of 
c-number function {fi(r)}\ 

«r)=Efli/<W (2-4) 

i 

with a,i and its hermitian conjugate a\ obeying the commutation relation 
[di, a]} = Sij, in accordance with (2.1). Take a normalized state vector | > 
which is also an eigenstate of all a\cii with 

a\di\ >= n t \ > . (2.5) 

For such a state, the expectation value of the Coulomb energy Ec ou i can be 
written as a sum of three terms: 

< \Hcoui\ >= E ex + E dir + E' dir (2.6) 

where 

E dir = d 3 rd 3 r'\ r - r'^ 1 < \J (r)\ X | J (r')| > (2.7) 

and 

^ = "E^/^V| r-rTVl/^r)! 2 !/^')! 2 - 

The last term is the subtraction, recognizing that in Wick's normal product 
each particle does not interact with itself. 

In the Schafroth solution, for the super phase at T = all particles are in 
the zero momentum state; therefore, on account of (2.2) the ensemble average 
of Jq is zero and so is the Coulomb energy. For the normal phase, take the 
magnetic field B = Bz with B uniform and pick its gauge field A = Bxy. At 
T = 0, let 

/i(r) = e^^x). (2.8) 
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Schafroth assumed pi = eBxi with Xi spaced at regular intervals A = 2ir/eBL, 
which approaches zero as L — > oo. This makes the boson density uniform and 
therefore E^ r = 0. In the same infinite volume limit, one can show readily 
that Q~ 1 E' dir — > 0. Since Schafroth omitted E ex , his energy consists only of 

E field = 



'field 

E-mech 



J <?r\ B\ (2.9) 



and 

r 1 

E d ia = Yl n iJ ^ T ^ l & ~ eA v( X )) 2 

= En l jd"r^(x-x l f. (2.11) 

The sum of (2.10) and (2.11) gives the usual cyclotron energy 

eB 

Emech + E dia = J2^i — . (2.12) 

Combining with (2.9), Schafroth derived the total Helmholtz free energy den- 
sity in the normal phase at zero temperature to be 

F n = - B 2 + ^- B (2.13) 
2 2m y J 

(Throughout the paper, we take e and B to be positive, since all energies are 
even in these parameters.) 

The derivation of (2.13) is, however, flawed by the omission of E ex . It turns 
out that for the above particle wave function (2.8), when Xi — Xj is << the 
cyclotron radius a = (e£>) 2, the coefficient of nirij in E ex is proportional to 
\xi — Xj\~ l . Hence Sl~ l E ex becomes 00 logarithmically as the spacing A — > 0. 

3. Corrected Normal State at High Density 

In this and the next section, we review the FLR analysis for the high density 
case, when p > r^ 3 where r& = Bohr radius = 47r/me 2 . 
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a. Strong field. We discuss first the case when B is >> (mp)*, so that 
the Coulomb correction to the magnetic energy (2.13) can be treated as a 
perturbation. To find the groundstate energy, we shall continue to assume 
(2.8) with pi = eBxi and X{ equally spaced at interval A, but keeping A ^ 0. 
Now as £1 — > oo, Vt~ l E' dir remains zero, but Q _1 Edi r in fact increases as A 2 
for A >> a, the cyclotron radius. The lowest value of E^r + E ex are both 
complicated in this range. The minimization can be done exactly, yielding 



A = 7raAo (3.1) 



where 



l-Q% = >ite-^'^ (3.2) 

The sum of E^r and E ex is found to be proportional to 1/B. Hence, (2.13) 
is replaced by 

with 

In = Iter + lex = 0.00567 + 0.00715 = 0.0128. (3.4) 



b. Weak field . Clearly (3.3) cannot be extended to B — > 0, as the last term 
would diverge. In its derivation the ipi(x) in (2.8) is taken to be the usual 
simple harmonic oscillator wave function determined by the magnetic field B 
only, without regard to Ecoul- This is valid only when B >> (mp)j. For 
much smaller B, we may consider a configuration in which the function ipi(x) 
is spread out flat over a width A — /, then drops to zero sinusoidally over a 
smaller width I on each side. Neighboring ipi overlap only in the strips of width 
/.Thus, it can be arranged that E |^(x)| 2 is uniform and hence 

i 

E dir = 0. (3.5) 

For sufficiently weak B, we find A >> the London length A^ = (e 2 p/my^. 
We must then drop the assumption that B is uniform; it is largest in the 
overlap region and drops to zero over the length A^ in either side. Let B 
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be the average of B over Q and a = (eB) 2 the corresponding "cyclotron 
radius". It is then found that 

A = =f , (3.6) 
I = \ L e 2 « \ L (3.7) 

where 

7T 5 \l/7 



is independent of B. The energies E mec h, E ex , E^ a are all proportional to B, 
and one obtains for the free energy density 

F n = \ B 2 + ^ rjB (3.9) 
2 m 

with 

. = ^»l, (3.10) 

much bigger than the Schafroth result, 
c. lntermediate field . Let 

B 1 = ep/m, B 2 = (mp)* (3.11) 

Between the above strong field case B >> B2 and the weak field case when 
B << r\B\, we have the regime when ip remains flat as in the previous section, 
but with \l » A > I. Hence, E^ r = as in (b), but B is uniform as in (a). 
One obtains an estimate 

F - = \ b2 + S b » /5b3/5 (3 - i2) 

where 

/16 \2/7 

B = Bl (3.13) 

The formulas (3.3-4), (3.9-10) and (3.12-13) are strictly upper bounds, 
which might be improved with better wave functions. We hope these to be 
good estimates. 
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4. Super State at High Density 



a. H Cl . The coherent length £ governing the disappearance of the normal 
phase outside a vortex is found from the Ginzburg-Landau (G.-L) equation[4] 
to be 

£ = (2A L /m)i (4.1) 

and hence (taking mXi >> 1) Xl >> £ so that we should have a type II 
superconductor, with the critical field for vortex penetration 

# Cl = ^i[ ln(A L /0 + 1.623 ] (4.2) 

in which the constant differs from that given by the G.-L. equation because of 
the long range Coulomb field. 

However, in Schafroth's solution because of (2.13), his normal phase would 
begin to exist at 

(H c )sch = \B 1 «H Cl , (4.3) 

above which [F n — BH)s c h would also be lower than that for the super phase, 
making it a Type I superconductor. But Schafroth's solution is invalid; (2.13) 
must be replaced by (3.9) at low B, giving 

H c = r)Bi (4.4) 

with rj >> \, as shown in (3.10). 

Next, we compare the above corrected H c with H Cl . Using (4.1), we write 
(4.2) as 

# Cl = (^lnC)£i (4.5) 

where 

lnC = | (§ln^ + 1.623). (4.6) 

Likewise, because of (3.8) and (3.10), the parameter 77 in (4.4) can also be 
expressed in terms of the same £: 

V = « C (4.7) 
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with the constant k, given by 



k = — 




2/7 



(4.8) 




(4.9) 



H C1 " 7 InC' 
Now C/l n C nas a minimum= e when ( = e. Thus, 




0.07 



> 1 



(4.10) 



and the system is indeed a Type II superconductor. 

b. Vortices . Once H = H Cl +, the vortices appear and soon become so 
numerous that their typical separation is of the order of Xl- This gives an 
average B of the order of B\ = ep/m. 

c. B > B\ . To increase B further, it is necessary to increase H on account of 
the interaction energy between vortices. The vortex separation distance is of 
the order of the cyclotron radius a = (eB)~*. In the regime £ << a « Xl 
(correspondingly, B 2 >> B » Bi), the vortices naturally form a lattice to 
minimize their interaction energy. An involved calculation gives the Helmholtz 
free energy density at T = to be 



d. B » £>2 - In this case of very strong magnetic field when the cyclotron 
radius and the separation distance are both much less than £ (but we assume 
that the system remains non-relativistic). The free energy density is dominated 
by the RHS of (2.13). However, there is still a super phase whose wave function 
is assumed to be given by the Abrikosov solution, giving by Eq.(8) of Ref.[5] 



F s = -B 2 + -£i£(ln -| + const) 
z 4 Jd 



(4.11) 



where B 2 = (rap) 2 and the constant is 




4.068 , square lattice 
4.048 , triangular lattice. 



(4.12) 
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and its Coulomb energy is calculated as a perturbation. The result for the 
Helmholtz energy density in the super phase is 

F s = -£ 2 + + (4.13) 

2 2m B 1 y ' 

with 

J 0.01405 , square lattice . . 

ls = \ 0.01099 , triangular lattice. ^ ' ' 

e. H C2 . The super phase regimes of Sections 4a-b, 4c, 4d correspond (with 
respect to the value of B or its average B) to those of the normal phase 
regimes discussed in 3b, 3c and 3a respectively. Since, for the triangular 
lattice, a comparison of (4.14) with (3.4) gives j s = 0.01099 < j n = 0.0128, 
we see that F s < F n for the same B » B2. Similarly, F s < F n for the same 
B in the regimes B2 » B » B\ and B « B\. From these results and 
that H = dF/clB is monotonic in B, one can readily deduce that the Legendre 
transform 

F = F — BE (4.15) 

satisfies F s < F n for the same H in all these regions. (See Figure 1) 
From this it seems possible that 

H C2 = 00; (4.16) 

the super phase may persist at high density for all values of the magnetic field. 

f. Remarks . In the problem discussed in Ref.[5], the Ginzburg-Landau function 
^ is an order parameter, whereas our fi(r) are single particle wave functions. 
Nevertheless, except for the constant in (4.2), the two problems have the same 
physics content at high p when H << B2. For higher field when H > B2, the 
Ginzburg-Landau ^ should vanish; however, this is not true in our problem. 
At T = 0, we place all the particles in the coherent state, making the charge 
density to vary greatly within a unit lattice cell. Our result (4.14) favoring 
a triangular lattice is unrelated to that of [5], because the ratio parameter 
< |\I/| 4 > / < |^| 2 > 2 in [5] does not appear in our problem. The lattice 
dependence in our problem is electrostatic in origin. 
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5. Low Density at Zero Field 

a. Normal Phase . At very low density and with zero magnetic field, E' dir of 
(2.7) becomes important. The lowest energy is now achieved by placing the 
individual charges in separate cells forming a lattice, with little or no overlap. 
Hence E ex can be disregarded, and a trial wave function leads in the limit 

p -> to 

e 2 

N-\E C oui + E mech ) = -o„— (5.1) 

where 

{4tt/3)R 3 = p' 1 and a n ^ 0.9 (5.2) 

very closely. The above formula (5.1) is valid for p « r b 3 , with r& the Bohr 
radius. 

b. Super Phase. In the same limit, the super phase energy also becomes 
negative, as shown by a Bogolubov-type transformation[6-8]. This leads to 

e 2 

N-\E Coul + E mech ) = -a s ^ ( 5 - 3 ) 

with 

0.316 <a s < 0.558. (5.4) 

Thus, a s < a n and the normal phase holds at p « r b 3 . 

c. Critical Density. As p increases, (5.1-2) serves only as a lower bound; i.e., 

N-\E Co ui + E mech ) > "(0.9)^. (5.5) 

For the normal phase, when p approaches r b 3 , the single particle wave function 
leading to (5.1-2) can no longer fit without overlap. We confine each particle 
within a cube, give it a r 1 sinqr wave function as a trial function, just avoid- 
ing overlap so that E ex = 0. With approximation neglecting the distinction 
between sphere and cube, we find 

where 

K n w 0.76 . (5.7) 
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Equating the above N^ 1 (Ec ou i + E mec h) for the normal phase with the corre- 
sponding expression (5.3) for the super phase, we find the critical density p c 
given by 

rl Pc = % (K n - a s f. (5.8) 

The system is in the normal state when p < p c , and in the super state when 
p > p c . (Eq.(4.12) in the FLR paper is equivalent to (5.8), but without the 
subtraction of K n by a s .) 

6. Further Improvement 

Although the FLR paper (66 pages in the Annals of Phys.) is quite lengthy, 
several important questions remain open. 

a. The energies in above sections 3 and 4 are all upper bounds obtained 
from trial functions. Perhaps a better trial function, like changing slabs into 
cylinders, might lower these bounds and put into questions some of the FLR 
conclusions. Also a numerical calculation exploring the transition regions 
would be valuable in case there are surprises, particularly when B ~ B^. 
In this connection we note that, e.g., in (4.9) the relevant factor in C/l n C 
is (mA^) 2 / 7 / ln(raAx), which becomes large when tuXl — > oo; yet, it is < 1 
when tuXl = 100 and only near but still less than 2 when ijiXl is 2000. 

b. The calculation of the above (5.6-7), i.e., Section 4.2 in the FLR paper, 
can be improved in several ways. First, consider the integral \ \ JqVSt, with 
V the potential due to Jo. Because the spatial integral of Jo is zero, and since 
each particle does not interact with itself we have 

N-\E dir + E' dir ) = e J ^ 2 (r)F(r)d 3 r (6.1) 

where ipo(r) is located inside a sphere, centered at zero, and V(r) is due 
to all of Jo except the term due to ipQ. Second, there is no need to ignore 
the distinction between sphere and cube. Using theorems from electrostatics, 
one can reduce (6.1) to the solution of a Madelung problem with like charges 
at lattice points and a background charge filling aJl space, plus a correction 
\e 2 p J ?/>o(r)r 2 <i 3 r. This correction can be combined with E mec h to optimize ipo, 
and the Madelung problem can be done by known methods. Third, the energy 
can probably be reduced by placing the centers of the particle wave functions 
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on a body-centered cubic lattice, as the cell available to each particle would 
then be more nearly spherical than a cube. 

c. Both FLR and the present paper have left open the question of what happens 
at low density and high field. It would be surprising if the boundary between 
super and normal phases were independent of the magnetic field strength. In 
the H versus p phase diagram at T = and high H, does the boundary 
between normal and super phases bend towards lower p, or towards higher? 

7. Comment 

The two most striking results in our paper are H Cl < H c , making the 
superconductor Type II instead of Type I, and that H C2 might be infinite. An 
improvement in the weak field normal trial function (the above Section 3b) 
might invalidate the first conclusion by lowering rj in (3.9). An improvement in 
the strong field normal trial function (Section 3a) could invalidate the second 
conclusion by lowering j n in (3.3). 

The field of condensed matter physics has received from its very beginning 
many deep and beautiful contributions from Russian physicists and masters L. 
D. Landau, V. L. Ginzburg, N. N. Bogolubov, A. A. Abrikosov, A. M. Polyakov 
and others. It is our privilege to add this small piece to honor this great and 
strong tradition and to celebrate the 90th birthday of V. L. Ginzburg. 
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Figure 1: Graphical construction of F = F — HB. At any point P on the curve F(B), 
the intercept of its tangent with ordinate gives F, since the tangent has a slope H. The 
subscript P denotes the values of B and F at P. 
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